Measurement of ep epn^n : experimental procedures and baryon resonance analysis 
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The cross section for the reaction ep e'pTT'^7v~ was measured in the resonance region for 
1.4<W<2.1 GeV and 0.5< <1.5 GeVVc^ 

using the CLAS detector at Jefferson Laboratory. 
The data show resonant structures not visible in previous experiments. The comparison of our data 
to a phenomenological prediction using available information on A'^* and A states shows an evident 
discrepancy. A better description of the data is obtained either by a sizeable change of the properties 
of the Pi3(1720) resonance or by introducing a new baryon state, not reported in published analyses. 

PACS numbers: 13.60.Le, 13.40.Gp, 14.20.Gk 



Electromagnetic excitation of nucleon resonances is 
sensitive to the spin and spatial structure of the transi- 
tion, which in turn is connected to fundamental proper- 
ties of baryon structure, like spin-flavor symmetries, con- 
finement, and effective degrees of freedom. In the mass 
region above 1.6 GeV, many overlapping baryon states 
are present, and some of them are not well known; mea- 
surement of the transition form factors of these states 
is important for our understanding of the internal dy- 
namics of baryons. Many of these high-mass excited 
states tend to decouple from the single-meson channels 
and to decay predominantly into multi-pion channels, 
such as Att or Np, leading to Nnn final states 0. 
Moreover, quark models with approximate (or "broken") 
SU(6)(g)0(3) symmetry 0, Q predict more states than 
have been found experimentally; QCD mixing effects 
could decouple these unobserved states from the pion- 
nucleon channel 21 while strongly coupling them to two- 
pion channels 0,0, . These states would therefore not 
be observable in reactions with ttN in the initial or final 
state. Other models, with different symmetry properties 
and a reduced number of degrees of freedom, as e.g. in 
ref. 0, predict fewer states. Experimental searches for at 
least some of the "missing" states predicted by the sym- 
metric quark models, which are not predicted by models 
using alternative symmetries, are crucial in discriminat- 
ing between these models. Electromagnetic amplitudes 
for some missing states are predicted to be sizeable ^ as 
well. Therefore, exclusive double-pion electroproduction 
is a fundamental tool in measuring poorly known states 
and possibly observing new ones. 

In this paper we report a measurement of the ep 
e'pTT~^Tr~ reaction studied with the CEBAF Large Ac- 
ceptance Spectrometer (CLAS) at Jefferson Lab. CLAS 
consists of a six-coil superconducting magnet producing 
an approximately toroidal magnetic field, allowing de- 
tection of electrons and hadrons and full 4-momentum 
recontruction. Three sets of drift chambers allow the de- 



termination of the momenta of the charged particles with 
polar angles from 10° to 140°. A complete coverage of 
scintillators allows the discrimination of particles by a 
time-of-fiight technique described in ref. [3|. The mag- 
netic field was set to bent positive particles outwards, 
away from the primary beam. We analyzed data taken 
in the so-called 'elc' running period, corresponding to 
about two months of data taking in the spring of 1999. 
Beam currents of a few uA were delivered to Hall B on 
a liquid- hydrogen target, corresponding to luminosities 
up to 4 X 10'^'^ cm~^s~^. The beam energies selected for 
this analysis were: 2.567 GeV at a torus current of 1500 
A ('low' field), with the goal of obtaining two bins, 
between 0.5 and 0.8 GeV^/c^ and between 0.8 and 1.1 
GeV^/c^, with W up to 1.9 GeV; 4.247 GeV at a field 
of 2250 A ('low' field), to obtain between 1.1 and 
1.5 GeV^/c^ , with ly up to 2.1 GeV. Important fea- 
tures of the CLAS 8j are its large kinematic coverage for 
multi-charged-particle final states and its good momen- 
tum resolution {Ap/p ~1%). Using an inclusive electron 
trigger, many exclusive hadronic final states were mea- 
sured simultaneously. Scattered electrons were identi- 
fied through cuts on the calorimeter energy loss and the 
Cerenkov photo-electron distribution. Different channels 
were separated through particle identification using time- 
of-flight information and other kinematic cuts. 

The hardware trigger in CLAS was based on a coin- 
cidence between Cherenkov counter and electromagnetic 
forward calorimeter. The discriminator threshold in the 
Cherenkov was put at a signal level of less than a single 
photoelectron produced on the photocathode, to avoid 
losing good events while rejecting most of the pion back- 
ground. The energy threshold in the calorimeter was set 
to cut off events as close as possible to the kinematic 
edges of the W and domain covered by the measure- 
ment, but sufficiently high to keep the low energy con- 
tamination due to hadronically interacting particles at a 
reasonable limit. A cut on the energy deposited in the 
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calorimeter was applied to eliminate part of the remain- 
ing pion background, mostly due to knock-off electrons 
giving a signal in the Cherenkov counter. We studied 
the effect of this cut on the electron detection efficiency, 
deriving a correction that was applied to the data. The 
Cherenkov efficiency was also studied by means of the 
shape of the photoelectron distribution in each module. 
This way, we evaluated the effect of a cut on the photo- 
electron number on the electron detection efficiency on 
each single Cherenkov module, as well as the module-to- 
module uniformity, deriving the appropriate corrections 
to be applied to the data. 

Another cut was applied in the (3 versus momentum 
plot, to select pions and protons in the reaction ana- 
lyzed. All Time of Flight (ToF) scintillator paddles were 
carefully inspected to make sure that the cut applied was 
appropriate for all of them. To identify the e'pTr+Tr" fi- 
nal state, we used the missing-mass technique, requiring 
detection in CLAS of at least ep7r+ . The good resolution 
allowed selection of the exclusive final state, ep-K^ir^ . 
After applying all cuts, our data sample included about 
2 X 10^ two-pion events. Fig. ^ shows the missing mass 
distribution obtained from CLAS, together with the ap- 
plied cut. The tail at higher mass is due to radiative 
effects and to multiple pion production. 




M/CGeVVc^) 

FIG. 1: Missing mass for detection of epix^ , at W^=1.6-1.7 
GeV and at =0.5-0.8 (GeV/c)^. Tlie tail at higher mass is 
due to radiative effects and to multiple pion production. The 
hatched area represents the adopted cut. 

In order to check the stability of CLAS in the de- 
tection of different reactions, we defined a set of his- 
tograms representing various reaction yields, i.e. elec- 
tron inclusive, electron-proton inclusive, electron elastic 
(W cut), electron-proton elastic (W plus 0-(j) correlation 
cut), electron-proton from the A (W cut) and finally 
electron-proton- 7r+ with a cut on the missing tt", each 
single yield being normalised to the Faraday cup charge 
obtained from the so called "live-gated" signal, where the 
signal from the Faraday cup is only integrated during the 
live-time: this way, the charge is already corrected for 
the data acquisition dead-time. Such normalised yields 
were calculated, using the PID procedures described in 
the previous subsection, for each data file inside a run 



(a run being typically a data taking over a period of an 
hour), therefore providing a very accurate monitoring of 
stability even inside a single run. "Good files" were se- 
lected requiring that the normalised yield for a single file 
should not deviate more than a few percent from the av- 
erage. Subsequent analysis was performed only on the 
good files. 

To obtain the cross section from the raw data, it is 
necessary to correct for detector non-uniformity, which 
has origin both in the geometry and in the response of 
the equipment to different particles. The geometrical and 
kinematic non-uniformity can be very well described by 
means of fiducial cuts that describe regions of the detec- 
tor where the response of the various subsystems is well 
known. Fiducial cuts will eliminate dead regions like the 
torus coils in the first place. The detector response inside 
the fiducial regions was simulated using a GEANT-based 
representation of the detector containing a detailed de- 
scription of particle interaction with the various subsys- 
tems. To evaluate all detector corrections, we divided 
the particle yield into kinematic bins defined through a 
complete set of independent kinematic variables of the 
hadronic state, plus W and . This corresponds to bin- 
ning all four-momenta of the particles involved (with the 
exception of the electron, for which no binning is applied 
to <j) azimuthal angle for symmetry reasons). This way, 
the folding of the detector response to the cross sections 
is done in narrow kinematic regions, where the cross sec- 
tion variation will be limited. The dependence of the 
results on the assumed particle distributions will there- 
fore be much reduced. However, events still have to be 
generated according to a realistic Monte Carlo, in order 
to minimize the model dependence of acceptance and ef- 
ficiency. 

The event generator used for simulations in this exper- 
iment contains several electroproduction cross sections, 
including single, double and triple pion electroproduc- 
tion. The code relies on cross section tables that describe 
measured total and differential cross sections from the lit- 
erature, scaled by a virtual photon flux and a dipole form 
factor to provide a reasonable fall off with . Therefore 
our code gives a realistic description of cross sections and 
their relative weights, as well as backgrounds generated 
from competing channels. Radiative effects were also in- 
cluded in the simulation. As an example of how the 
Monte Carlo reproduces the main features seen in the 
data, we reproduce here the comparison of some kine- 
matic distributions from the data at 2.567 GeV beam 
energy, 1500 A field, with the corresponding simulation. 
As specific example we report the case of the W bin 1.6- 
1.625 GeV, and between 0.5 and 0.8 GeVVc^. Fig- 
ure |21 shows the invariant mass distribution for the pion- 
pion pair. Figure |31 shows the invariant mass distribution 
for the proton- 7r+ pair. Figure^ shows the CM angle of 
the p — TT^ pair (which would correspond to the CM angle 
of a A° in the specific case that a A° is produced). 

A particularly important issue regards the percentage 
of events lost in the binning process due to bins with zero 
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FIG. 2: Simulated (top) and measured (bottom) invariant 
mass for the pion-pion pair, when detecting epir^ , at 2.567 
GeV beam energy, 1500 A field, W between 1.6 and 1.625 
GeV and between 0.5 and 0.8 GeVVc^- 




100 120 140 160 180 

t5^"p,- (deg) 



TlietaDO for Data: 2pii 



t5-"p,_ (deg) 



FIG. 4; Simulated (top) and measured (bottom) angular dis- 
tribution of the p — n~ pair in the CM system, when detecting 
ep-K^ , at 2.567 GeV beam energy, 1500 A field, W between 
1.6 and 1.625 GeV and between 0.5 and 0.8 GeVVc^ 




FIG. 3: Simulated (top) and measured (bottom) invariant 
mass for the proton- 7r+ pair, when detecting epTr"*", at 2.567 
GeV beam energy, 1500 A field, W between 1.6 and 1.625 
GeV and between 0.5 and 0.8 GeVVc^- 



acceptance: of course, for such bins it is not possible to 
correct the data and obtain a cross section; being the 
kinematic variable space multidimensional in the case of 
double pion production, extrapolating the cross sections 
from neighboring bins with non-zero acceptance to bins 
vifith zero acceptance can be unreliable. Therefore, we 
checked carefully, using our realistic Monte Carlo, the 
percentage of cross section lost because of bins with zero 
acceptance. It turned out that with the adopted binning, 
the percentage of unmeasured cross section in bins with 
vanishing acceptance was below 10 % and typically of the 
order of a few percent. Actually, the fact that a cell has 
zero acceptance or efhciency may be just connected to 
insufficient statistics in the simulation, as many cells end 
up with very few generated events; nevertheless, for those 
cells it is not possible to perform a correction to the data. 



To get the final cross sections, we therefore performed an 
extrapolation to the empty bins using the MonteCarlo 
distributions as an estimate of the missing cross section. 
This extrapolation is typically only a few percent and we 
quoted as systematic error a quantity equal to one half of 
the extrapolation, assuming that the extrapolated cross 
section may be wrong by plus or minus 50 %, due to the 
assumptions in the MonteCarlo. 

The range of invariant hadronic center-of-mass (CM) 
energy W (in 25 MeV bins) was 1.4-1.9 GeV for the 
first two bins in the invariant momentum transfer Q^, 
0.5-0.8 (GeV/c)2 and 0.8-1.1 (GeV/c)^, and 1.4-2.1 GeV 
for the highest bin, 1.1-1.5 (GeV/c)^. Following the 
procedures schematically described above, data were cor- 
rected for acceptance, reconstruction efficiency, radiative 
effects, and empty target counts. They were further 
binned in the following set of hadronic CM variables: in- 
variant mass of the pTt"*" pair (10 bins), invariant mass 
of the 7r+7r~ pair (10 bins), tt~ polar angle 9 (10 bins), 
aziniuthal angle (f> (5 bins), and rotation freedom ip of 
the pT:~^ pair with respect to the hadronic plane (5 bins). 
The fully differential cross section is of the form: 



da 



dWdQ^dMp^+dM^+^-d COS 0„- d(j)^- d^p^4 



d<T„ 



(iMp^+dM^+^- d cos 6*^- dcf).,:- dijjp 



dffy 



47r E^M^ 



1 WiW"^ - M2) 
(l-e)g2 



(1) 



(2) 



where r„ is the virtual photon flux, is the virtual 
photon cross section, a is the fine structure constant, E 
is the electron beam energy, Mp is the proton mass, and 
e is the virtual photon transverse polarization 9J. 

Since existing theoretical models 11] are limited to 
W <1.6 GeV, we have employed a phenomenological 
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CM (deg) 

from CLAS (from top 



to bottom) at W=1.5-1.525 GeV and for the three mentioned 
intervals (left to right). The error bars include statisti- 
cal errors only. The curves represent our step (A) reference 
calculations. 
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FIG. 7: 



CM (deg) 

J"- , and ^ from CLAS (from top 

to bottom) at W^=l. 625-1. 65 GeV and for the three mentioned 
intervals (left to right). The error bars include statistical 
errors only. The curves represent our step (A) reference cal- 
culations. 
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and 



CM (deg) 

^ from CLAS (from top 



to bottom) at W=1.575-1.6 GeV and for the three mentioned 
intervals (left to right). The error bars include statisti- 
cal errors only. The curves represent our step (A) reference 
calculations. 
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FIG. 8: 



and 



CM (deg) 

from CLAS (from top 



to bottom) at VK=1.775-L8 GeV and for the three mentioned 
intervals (left to right). The error bars include statisti- 
cal errors only. The curves represent our step (A) reference 
calculations. 



calculation [l^ for a first interpretation of the data. 
This model describes the reaction ^^p — > ptt^tt" in the 
kinematic range of interest as a sum of amplitudes for 
7t,p — > Att pTT+TT^ and ^^p — > p'^p — > ptt+tt", while 
all other possible mechanisms are parameterized as phase 
space. A detailed treatment was developed for the non- 
resonant contributions to Att, while for pp production 
they were described through a diffractive ansatz. For the 
resonant part, a total of 12 states, classified as 4* Q], with 
sizeable Att and/or pp decays, were included based on a 
Breit-Wigner ansatz. A few model parameters in non- 
resonant production were fitted to CLAS data at high 
W , where the non-resonant part creates a forward peak- 
ing in the angular distributions, and kept fixed in the 



subsequent analysis. The phase between resonant and 
non-resonant Atf mechanisms was fitted to the CLAS 
data as well. To simplify the fits, we reduced eqn. (1) to 
three single-differential cross sections, — , -rrr^ — , 

and -j-^^^fg — , by integrating over the other hadronic vari- 
ables. These three 1-D distributions were then fitted si- 
multaneously. Here -r-^ — and . — are both con- 

nected with the dominant A^+tt" production reaction, 
while -jf-f^ — is connected with pp^ production. For 

each W and bin, a total of 26 data points from the 
three single-differential cross sections were used in our 
fits. The two edge points in both the pTr"*" and 7r"'"7r~ 
mass distributions were excluded as the model did not 
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FIG. 9: 



CM (deg) 

and , '^"^ from CLAS (rrom top 

to bottom) at W^=l. 825-1. 85 GeV and for the three mentioned 
intervals (left to right). The error bars include statistical 
errors only. The curves represent our step (A) reference cal- 
culations. 
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FIG. 10: Left: Total cross section for -fvP — * ptt^tt 
as a function of W . Data from CLAS are shown at 
Q2=0.5-0.8 (GeV/c)2 (fuU points), ^^^0.8-1.1 {GeV/cf 
(open squares), and Q^=l.l-1.5 (GeV/c)'^ (open triangles). 
Error bars are statistical only, while the bottom band shows 
the sytematic error for the lowest bin. The curves rep- 
resent our step (A) reference calculations. Right: ^^^f"^ — 

from CLAS at Q2=0.8-1.1 (GeV/c)^ and W=1.7-1.725''GeV 
(statistical error bars only). The curves represent our step 
(A) reference calculations, extrapolated to the edge points. 
The dashed line includes all resonances, the dot-dashed line 
includes only the non-resonant part, and the solid line is the 
full calculation. 



take into account the kinematic smearing in the M^+t^- 
versus Mp^+ plot caused by the W bin width. 

The physics analysis included the following steps: (A) 
We produced reference curves using the available infor- 
mation on the iV* and A resonances in 1.2-2 GeV mass 
range. Discrepancies between the CLAS data and our 
calculation were observed, which led to the subsequent 
steps B and C. (B) Data around W~\.l GeV were fit- 
ted using the known resonances in the PDG but allowing 



the resonance parameters to vary in a number of ways. 
The best fit, corresponding to a prominent P13 partial 
wave, could be attributed to the PDG Pi3(1720) reso- 
nance, but with parameters significantly modified from 
the PDG values. (C) As an alternative description, we 
introduced a new baryon state around 1.7 GeV. In what 
follows we describe each of the steps above in more detail. 

Step (A) - To produce our reference curves, the evo- 
lution of the A112 and ^3/2 electromagnetic couplings for 
the states was taken either from parameterizations of ex- 
isting data Il3|. or from Single Quark Transition Model 
(SQTM) fits 'IS'l where no data were available. For the 
Pii(1440) (Roper), given the scarce available data, the 
amplitude was taken from a Non-Relativistic Quark 
Model (NRQM) jlj]. Partial LS decay widths were taken 
from a previous analysis of hadronic data and renor- 
malized to the total widths from Ref. Results from 
step (A) are reported in Figs.[Slto|51for the specific one- 
dimensional differential cross sections analysed and for a 
sample of W bins. In Fig. ^| (left), we report the total 
cross section data and the corresponding curves from step 
(A). As demonstrated by the plots, our calculation, even 
without performing any adjustment of the resonance pa- 
rameters, is able to give a good accounting of the main 
features of the data in a wide W and region. The 
total cross section strength for W < 1.65 GeV (except 
for the region close to threshold), and ior W > 1.8 GeV 
is well reproduced. However, a strong discrepancy is ev- 
ident at W around 1.7 GeV. Moreover, at this energy 
the reference curve exhibits a strong peak in the tt+tt^ 
invariant mass (Fig. [TUl right), connected to sizeable p 
meson production. This contribution was traced back 
to the 70-91% branching ratio of the Pi3(1720) into this 
channel Q, 111 111. 

Step (B) - We then considered whether the observed 
discrepancy around 1.7 GeV could be accomodated by 
varying the electromagnetic excitation of one or more of 
the PDG states. In our code, interference effects between 
continuum and resonances like those reported in Ref. , 
are taken into account automatically inside the model 
frame Therefore, our investigation at this stage was 
including the possibility of accounting for the 1.7 GeV 
structure via interference effects, although the peaking 
of such an interference pattern at the same W for all 
bins would be rather surprising. All fit values were 
calculated from the 8 W bins between 1.64 and 1.81 GeV 
and from the 3 bins (624 data points). The number 
of free parameters ranged from 11 to 32, depending on 
the fit, corresponding to i'—6l3 to 592 degrees of free- 
dom. Assuming the resonance properties given by the 
PDG, the bump at about W—1.7 GeV cannot be due to 
the 1515(1675), Fi5(1680), or £133(1700) states; the first 
because its well known position cannot match the peak; 
the second because of its well known position and pho- 
tocouplings the third due to its large width (^300 
MeV). The remaining possibilities from the PDG are the 
£)i3(1700), the Pi3(1720), and the Pii(1710) (the latter 
was not included in step (A)), as there are no data avail- 
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TABLE I: PDG Pi3(1720) parameters from fit (B) and new 
state parameters from fit (C). Errors are statistical. 





M (MeV) 


F (MeV) 


^ (%) 


(%) 


PDG Pi3 (B) 


1725±20 


114±19 


63±12 


19±9 


PDG fl] 


1650-1750 


100-200 


N/A 


70-85 


new Pi3 (C) 


1720±20 


88±17 


41±13 


17±10 



able on the dependence of A1/2 or lUf- If 

configuration mixing occurs, the Z3i3(1700) cannot be ex- 
cited in the SQTM from proton targets, while the SQTM 
prediction for the P13 (1720) relies on ad hoc assumptions. 
According to the literature hadronic couplings 

of the 1313(1700) and the total width of the Pii(1710) are 
poorly known, while the Pi3(1720) hadronic parameters 
should be better established. Several other partial waves 
were investigated in step (C). 



Pii(1710). In both fits (B2) and (B3), we also varied the 
hadronic parameters and the position of the Di3(1700) 
over a range consistent with their large uncertainties from 
Ref. m. Fits (Bl) and (B3) gave a poor description of 
the data, with x^/i'—5.2 and 4.3, respectively. The best 
fit (x'^/iy =3.4) was obtained in (B2) (Fig. [ill). How- 
ever, the resulting values for the branching fractions of 
the P 13 (1720) were significantly different from previous 
analyses reported in the literature and well outside the re- 
ported errors dllllll. Starting from (B3), we then per- 
formed a final fit, (B4), for which the Pi3(1720) hadronic 
couplings were fixed from the literature, and varying the 
photocouphngs of all three candidate states, Di3(1700), 
Pi3(1720), and Pii(1710), by 100% {a). No better solu- 
tion was found, the x^/^ being 4.3 (Fig. [TT)l. In Fig. IT^ 
we report the final comparison of fits (B2) and (B4) with 
the total cross section data. Table shows our results 





FIG. 11 



to bottom) at 11^=1.7-1.725 GeV and for the three mentioned 
intervals (left to right). The error bars include statisti- 
cal errors only. Curves (see text) correspond to the fits (B2) 
(red), and (B4) (blue), and are extrapolated to the mass dis- 
tributions edge points. 

To improve our reference curves before fitting the 
bump at around 1.7 GeV, the following steps were car- 
ried out: the Pii(1440) strength was fitted to our low W 
data; the 1)15(1675) and the Z)i3(1700) photocouphngs 
(which vanish in the SQTM) were replaced by NRQM 
values from Ref. an empirically established Ai/2,3/2 
SQTM fitting uncertainty or NRQM uncertainty of 10 
or 20% (ct) was applied to all N* states; the hadronic 
parameters were allowed to vary for the 1)13(1700) ac- 
cording to Ref. ^^; and finally, the curves providing 
the best x^/z^ were selected as the starting points. First 
we performed three fits, (Bl), (B2), and (B3), where the 
photocouphngs of only one resonance at a time were var- 
ied. In (Bl), we varied A1/2, A3/2, hadronic couplings, 
and position of the Di3(1700) in a wide range. In (B2), 
the same was done for the Pi3(1720), and in (B3) for the 



1.6 1.8 

W (GeV) 

FIG. 12: Left: Total cross section for 'y^p — > pn'^Tv' as a 
function of W from CLAS at the 3 mentioned intervals 
(see fig.l). The error bars are statistical only. The curves 
(see text) correspond to the fits (B2) (red), and (B4) (blue). 
Right: subdivision of the fitted cross section (B2) for Q^=0.5- 
from CLAS (from top g.S (GeV/c)^ into resonant A++7r- (black solid), continuum 
A^"'"7r~ (black dashed), resonant p^p (magenta solid), and 



100 

CM (deg) 



continuum p^p (magenta dashed). Notice the different verti- 
cal scales. 



(first row) with statistical uncertainties, in comparison 
with the PDG values (second row). Our fits were not 
providing an unambigous separation of ^1/2, ^3/21 and 
the longitudinal 5'i/2, so we report as result the total 

photocoupling strength. 



"^3/2 



Sf/^. The re- 



sulting value for the Pi3(1720) fit is reported in the first 
three rows of Table |n] The errors reflect the statistical 
uncertainties in the data and the correlations among the 
different resonances. 

As discussed above, fitting the data around 1.7 GeV 
with established baryon states leads either to a poor fit or 
to a drastic change in resonance parameters with respect 
to published results. In the framework of our analysis, 
there is no way to assess the reliability of the previously 
determined hadronic parameters of the PDG Pi3(1720). 
The resonant content of the reaction ttN — > ttttN, which 
is used to obtain the hadronic parameters, may be dif- 
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TABLE II: PDG Pi3(1720) total photocoupling from fit (B2) 
and new state total photocoupling from fit (C). Errors are 
statistical. 



step 


(GeV/c)2 


\l ^1/2 + ^3/2 + '^1/2 

(10-VVGeV) 


B2 


0.65 


83±5 


B2 


0.95 


63±8 


B2 


1.30 


45±27 


C 


0.65 


76±9 


C 


0.95 


54±7 


C 


1.30 


41±18 



ferent from that of reactions initiated by an electromag- 
netic probe. In particular, the Pi3(1720) state seen in 
TT-ZV — > TTTriV may not be excited in electroproduction, 
while some other state that decouples from -kN may be 
excited electromagnetically. This possibility is studied in 
the next step. 

Step (C) - We investigated whether our data could be 
fitted by including another baryon state, while keeping 
the hadronic parameters of the Pi3(1720) as in Refs. 
P, 10. The quantum numbers 5/i,P7i,P73,D/3,Z?75, 
F/5,F/7 were tested on an equal footing, where //2 is the 
isospin, undetermined in our measurement. We then si- 
multaneously varied the photocouplings and the hadronic 
parameters of the new state and the Di3(1700). The to- 
tal decay width of the new state was varied in the range of 
40-600 MeV, while its position was varied from 1.68-1.76 
GeV. The best fit (x^/i^=3.3) was obtained with a P/3 
state, while keeping the Pi3(1720) hadronic parameters 
at published values. Other partial waves gave a. /v > 
4.2. Curves obtained from the best fit were nearly identi- 
cal with the red sohd lines in Figs.llllandll2l In order to 
avoid the unobserved p production peak (Fig. ^| right), 
the photocouplings of the PDG Pi3(1720) had to be re- 
duced by about a factor of two with respect to the SQTM 
prediction, making its contribution very small. Instead, 
in this fit the main contribution to the bump came from 
the new state. Resonance parameters and total photo- 
coupling value obtained from the assumed new state are 
reported in Table Q (last row) and UTI flast 3 rows), re- 
spectively. 

At this point, one could wonder whether a similar dis- 
crepancy between the calculation and the data is present 
in the real photon case, and whether introduction of a 
new state is compatible with the existing real photon 
data as well. To this purpose, we reproduce in Fig. El 
the old data from Ref. [Ig , together with curves obtained 
from our model. The red solid curve in the graph repre- 
sents the overall contribution from the non-resonant pro- 
cesses, while the blue dashed line corresponds the res- 
onant contribution only, without introducing any new 
state. The black dashed line shows the result of our 
model when including both non-resonant and resonant 
processes, without new states. Finally, the black solid 



curve shows our fit of the real photon data when in- 
troducing a new state, and the blue solid curve shows 
the corresponding modification in the pure resonant ex- 
citation curve. The photocouplings of the new |^(1720) 
state were determined extrapolating the values obtained 
from the CLAS experiment with virtual photons. An 
additional adjustment, within 10 %, was applied to re- 
produce the invariant mass distributions of the pairs pn'^ 
and TT+TT", measured at W=1.7 GeV It is clear how, 
due to the stronger dominance of the non-resonant mech- 
anisms at the photon point, as well as due to the partic- 
ular interference effects, a new state can be accomodated 
into the picture. In fact, the agreement between data 
and our calculation is even improved by the introduction 
of the new state. 
100 I 1 




1.4 1.5 1.6 1.7 1.8 1.9 2 

W (GeV) 

FIG. 13: Total cross section for 7p —> pn^-n~ as a function 
of W from Ref. (3j together with curves obtained from our 
model. Red solid line: overall contribution from the non- 
resonant processes; blue dashed line: resonant contribution 
only, without introducing any new state; black dashed line: 
full calculation, without new states; black solid curve: full cal- 
culation, including the new state; blue solid curve: resonant 
contribution only, including the new state. 

A second P13 state was indeed predicted in Ref. Q, 
with a mass of 1870 MeV, and in Ref. with a 

mass of 1816 MeV. The presence of a new three-quark 
state with the same quantum numbers as the conven- 
tional Pi3(1720) in the same mass range would likely lead 
to strong mixing. However, as mentioned above, a dif- 
ferent isospin and/or partial wave cannot be excluded. 
The new state may also have a different internal struc- 
ture, such as a hybrid baryon with excited glue compo- 
nents. Such a P13 hybrid state is predicted in the flux 
tube model 20]. Finally, in Ref. '21] the existence of a 
P13 pentaquark (qqqqq) configuration with a mass in the 
range 1.76-1.78 GeV is predicted. Yet another possibility 
is that some resonance parameters established in previ- 
ous analyses may have much larger uncertainties than 
reported in the literature. In this case, outlined in our 
step (B), our analysis would establish new, more precise 
parameters for a known state, and invalidate previous 
results. 

In conclusion, in this paper we presented data on the 
ep — > e'pTT+TT^ reaction in a wide kinematic range, with 
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higher quality than any previous double pion produc- 
tion experiment. Our phenomenological calculations us- 
ing existing PDG parameters provided a poor agreement 
with the new data at ~ 1700 MeV. We explored 
two alternative interpretations of the data. If we dis- 
miss previously established hadronic parameters for the 
Pi3(1720) we can fit the data with a state having the 
same spin/parity/isospin but strongly different hadronic 
couplings from the PDG state. If, alternatively, we in- 
troduce a new state in addition to the PDG state with 
about the same mass, spin |, and positive parity, a good 
fit is obtained for a state having a rather narrow width, a 
strong Att coupling, and a small pN coupling, while keep- 
ing the PDG Pi3(1720) hadronic parameters at published 
values. In either case we determined the total photocou- 
pling at (5^ > 0. A simultaneous analysis of single and 
double-pion processes provides more constraints and may 



help discriminate better between alternative interpreta- 
tions of the observed resonance structure in the CLAS 
data. Such an effort is currently underway. 
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